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Dynamics of interacting quintessence 
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In this paper, we investigate coupled quintessence with scaling potential assuming specihc forms 
of the coupling as A namely, ap'm, P(><i> and (j{pm + p</>)i present phase space analysis for three 
different interacting models. We focus on the attractor solutions that can give rise to late time 
acceleration with UIde/Ddm of order unity in order to alleviate the coincidence problem. 


I. INTRODUCTION 

A large number of cosmological observations [H-Q reveal that our Universe is experiencing an accel¬ 
erated expansion at present and the transition from deceleration phase to acceleration phase took place 
in the recent past [3| . In the standard Einstein gravity, the late time cosmic acceleration is driven W 
an exotic energy component with huge negative pressure filling the Universe, known as ‘dorfc energy’’ [6|. 
One of the simplest candidate of dark energy (DE) is the cosmological constant (CC) A. However, it is 
plagued with difficult theoretical issues such as fine tuning and cosmic coincidence problem Q. This is 
important to explore whether dark energy is cosmological constant or it has dynamics. To this effect, 
a variety of dynamical dark energy models have been explored in the references [sl-l^. The models of 
unstable dark matter with a non zero cosmological constant can also mimic such dynamics [2lj. Al¬ 
ternatively, large scale modification of gravity has been used to obtain late time cosmic acceleration. 
At present late time cosmic acceleration is treated as an established phenomenon however its underly¬ 
ing cause is still unknown. Within the framework of Einstein gravity and modified theories of gravity, 
numerous models can explain the said phenomenon. 

Although ACDM model is consistent with present observations, yet there is no satisfactory argument 
for coincidence problem. Interaction of dark energy with dark matter is one novel approach that might 
address the mentioned problem. The interacting dark energy models have been recently proposed by 
several authors (22l - [^ . The interaction between dark energy and dark matter may enhance the dark 
matter, and also affect structure formation. The investigation of phase space analysis is the one conclusive 
test for dark energy models. Specifically, the attractor solutions are independent for a wide range of initial 
conditions. If the dark energy models have ^de/^dm of the order 1 and an accelerated scaling attractor 
solution, then the coincidence problem can be alleviated. The non-interacting quintessence psi . and 
quintom [s^ models show late time accelerated attractors, and possess UIde as 1, therefore, they do 
not provide an adequate solution for coincidence problem. In the literature, two forms of interactions 
have been discussed namely local and non- local. Local forms of interactions are directly proportional 
to energy density whereas non- local forms are directly proportional to Hubble parameter H and energy 
density p. In this paper we consider local forms of interactions proportional to energy density. Some of 
the local forms have been discussed in references [8ll - [^ . Note that some of the choices of interacting 
terms appeared implicitly in the literature [s^. There is also approach to discuss the interacting term 
without the assumption of a specific form of interacting term |36| . The plan of the work is organised as 
follows: In section m we establish the interacting quintessence cosmological framework and construct an 
autonomous dynamical system which is worthy for phase space investigation. In section IIIII we discuss 
phase space analysis and hnd stationary points and their stability for three interacting quintessence 
cosmological models. Our results are presented in section HVl 
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II. QUINTESSENCE COSMOLOGY 


We consider two components first one is canonical scalar field (quintessence) as a source of dark 
energy in spatially flat Universe, and second one is matter (Baryonic+DM). The total energy density of 
the Universe is conserved, and the individual components of energy density may not be conserved. Thus, 
we are considering following conservation equations of energy density as: 


Ptot + 3iL(l + Wtot)ptot — 0, 
p^ + 3i7(l + Wtj,)ptj, = -A, 

Pm T 3H[I T WrrCjPm — A^ (1) 


where ptot = P<f> + Pm, 1 + = 9^ / Ptf,, Wm = 0 is the equation of state of matter, A is the interaction 

strength and H is the Hubble parameter which is given as 




SttG 

^ Ptot 


( 2 ) 


The sign of A gives information about the direction of flow of energy between two components. There 
are 3 cases: 

Case I: If H > 0, in this case transfer of energy occurs from quintessence to dark matter. Consequently, 
quintessence losses self strength and gives dark matter. 

Case II: HA <0, under this condition dark matter losses its strength and there is energy transfer from 
dark matter to quintessence. 

Case III: If H = 0, under this condition quintessence do not interact with dark matter, and no energy 
transfer at all between two components considered in the literature. Therefore, we are not considering 
this case. 

Since, there is no fundamental theory of dark energy and dark matter interaction (interaction in dark 
sector) at present, therefore it is not possible to construct the functional form of interaction strength A 
from first principle. Different forms of interaction strength (linear and non-linear) have been considered 
by several authors [37l-l45j|. Motivated from the left hand side of the energy conservation equation (P) it 
is natural that A should be the function of Hubble parameter and energy density that is 


A = A{H, Pm,P<t>) 

Here we consider three specific forms of interaction strength heuristically as: 


( 3 ) 


A — cypjn, 

A = Pp^, 

A = a{pm + P4>) 


( 4 ) 

( 5 ) 

( 6 ) 


Since, Hubble parameter has the dimension of inverse of time, and inverse of time is sitting in the rate 
change of energy density, see equations (P)- ®, therefore we are not incorporating H separately in 
the functional form of A. The evolution equations in a flat Friedmann- Lemaitre- Robertson- Walker 
(FLRW) Universe can be written as: 


H — {pm + P4>) 


( 7 ) 


“21111 — —(pm+pV) 

where = SttG, p^ = + V{(j)) and = ^(j)^ — V{(j)). We introduce following dimensionless 

parameters 


^2 ^ 


= 


K^V 


A = -^ 
kV 


6i72 ’ * 3ip2’ 

to form an autonomous system of evolution equations m and © as: 


( 8 ) 


dX 

d/V 

dY 

d/V 


-3Y-t 





( 9 ) 
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where N = Ina. The total equation of state and field density parameter are given as: 


Wtot 


- 1 - 


2H 


3i?2 


= 


The condition for acceleration is Wtot < 


( 10 ) 


III. PHASE SPACE ANALYSIS: STATIONARY POINTS AND THEIR STABILITY 

In this section we shall use an autonomous system [equation ([H])], which is appropriate for obtaining 
stationary points and stability. The stationary points shall be obtained by equating the left hand side 
of equation to zero. Their stability will be confirmed from the sign of the corresponding eigenvalues, 
which will be obtained numerically. 


A. Interacting model I 


In this model we consider the following specific form of interaction strength 

A = apm ( 11 ) 

Using equation © for this phenomenological form of interaction we have 

H _ S 

Ip ~~2 


l-Y^ + X^{l-2a) 


1 — a 


( 12 ) 


By using +17^ = 1 where 17^, = 1 — — Y^. Thus, from equation ® with equation (TT^ we have 

following form of autonomous system 


dX 

dN 

dN 


= -3A + 


i/^AY^ + ^X 

'1-Y2 + X2(l-2a)' 

+2^ 

1 — a 


(13) 


\f^XXY + ^Y 

'1 - Y2 + X2(1 - 2a)' 

V 2^ + 2 

1 — a 


Using equations OT and (na, the total equation of state for this model can be written as 


Wtot 


-1 + 


'1 - Y^ +X^{1 - 2a)' 
1 — a 


= w^n^ 


(14) 


The critical points of autonomous system (USD could be obtained by setting = 0 and = 0 
simultaneously. Thus, we have following stationary points: 


(1) X = -1, Y = 0 

In this case, the corresponding eigenvalues are 

Pi = 6 + 3/(—1 + a) <0, fora <1, p 2 = 3 + \/3/2 A <0, for A <—'/E 

The eigenvalues of this point show the negativity for a < 1 and A < —i/E. Therefore, it is stable point. 

(2) X = 0, Y = 0 

In this case, the corresponding eigenvalues are 

Pi = 3/(2 — 2a) < 0, for a > 1, P 2 = —3 + 3/(2 — 2a) < 0, for 1 < a < 1/2 
The negativity of the eigenvalues represents the stability. This point is stable for a >1. 


(3) X = 1, Y = 0 

In this case, the corresponding eigenvalues are 
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(a) 



(c) 



FIG. 1: This figure shows the phase portrait, evolution of energy density and density parameter for model I for 
the stable fixed point 5 which is an attractive node. The panel (a) corresponds to a = —3.6 and A = 0.94, for 
these values of the parameters we obtain = 0.73, Wtot = —0.78, = —1.06 and an accelerating attractor 

solution. The panel (b) corresponds to a = —4.6 and A = 0.94, correspondingly we get = 0.60, Wtot = —0.82, 
= —1.36 and an accelerated attractor solution. In both the panels, black dots represent attractor stable 
point. The panels (c) and (d) have same values of the parameters as panel (a), and show scaling behaviours that 
provide an accelerated expansion. 


/Ti = 6 + 3/(—1 + a) <0, fora <1, ^2=3— •\/3/2 A <0, for A < 

This point is stable for a <1 and A < \/6. 




1 - 4 ^ 


“ V6’ " “ V " 6 

In this case, the corresponding eigenvalues are 
/Ti = (—6 + A^)/2 <0, for A < \/6, /X 2 = 3/(—1 + a) + A^ 

This point is stable under above given conditions. 


< 0 , 


for a <1, A < •\/3/(—1 + a) 


(5) X = 


A(l-a) ’ 


F = 




3a 


A(a-l) 

In this case, the corresponding eigenvalues are 


3(A^(a-l)^(2a-l) + .^A2(a-l)3(2a-l)(24+A2(7+2a)(a-l) ^ ^ 
Ml — 4A2(a-l)3 < 


for a < A > 0, 
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3(A^(a-l)^(2a-l)+^A2(a-l)3(2Q:-l)(24+A2(7+2a)(a-l) 

4A2(a-l)3 


for a < -7/2, A > 


This point has negative eigenvalues for a. < —712 and A > y^ 2/3. Therefore, it is stable point. We 
are interested in this point because it has both the parameters a and A. We evolve the autonomous 
system m numerically for the parameter values a = —3.6, A = 0.94, and a = —4.6, A = 0.94, and 
the obtained results are shown in figure [1] For the chosen parameters the stable point behaves as an 
attractive node which is confirmed by the panels (a) and (b) of figure [T] The lower panels of hgure 
[T] show the scaling behaviour that gives acceleration. In addition we also calculate basic cosmological 
observables Wtot, flfj, and obtained as —0.78, —1.06, 0.73 and —0.82, —1.36, 0.60 corresponding to 
a = —3.6, A = 0.94, and a = —4.6, A = 0.94, respectively. This point is summarized in table HI 


B. Interacting model II 


This model is specified by a coupling of the form 

A = /3p0 

For this coupling term we have 


H 

IP 

dX 

dN 

dY 

m 


3 

~2 

= -3X 


1-^2^X^(1-/3) 


1 + /3 


1 - 


X^l-p) 


1 + /3 




1 + /3 


Using equations (0 and (HU), the total equation of state for this model can be written as 


Wtot = -1 


1 -^ 2 ^^^( 1 -^) 


1 + /3 


= 


This interacting model has following stationary points: 


(15) 


(16) 

(17) 


(18) 


(1) X = 0, F = 0 

In this case, the corresponding eigenvalues are 
III = —3/2, /i 2 = 3/2 

This is unstable point because one of the eigenvalue is positive. 


( 2 )v = ^, y = o 

In this case, the corresponding eigenvalues are 

III = 3, ^2 = 3 - < 0 , for ^-3(l + /3)/2 A>3(-l+/3) 

It is also unstable point because its one eigenvalue is positive. 


^o^ (3 +aU(1+/3)-<5 _ V(1--S^)(6+9/AU-(1+/3^)A"+-5(1+/3-3/A2+3/3/A2) 

'-'^1 “ 2V^ A ’ “ 2^3 

In this case, the corresponding eigenvalues are 

111 = ^(3 (-5 + A^) + 3/3 (3 + A^) - 3 ,5 - V ( 2 I 6 - 63A2 + A^ + ^2 ^_24 y 2 ^ ^3 ^ 2^2 _ _ ^^ 2 ^ _ y4j_^ 

/3 (2A® + 72,5 - A‘^(18 + 6 ) + 3A2(30 + 7,5))))) 

112 = ^(3 (-5 + A2) + 3/3 (3 + A2) - 3 5 + V ( 2 I 6 - 63A2 + A^ + /32 (-24 + A^) (3 + A^)^ - 726 - 3A25 - A45+ 
/3 (2A® + 725 - A4(18 + 5) + 3A2(30 + 75))))), 


where 


5 = v^(l + /3) ((-3 + A2)2 + /3(3 + A2)2) 


(19) 
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(a) 



ln(l+z) 


(b) 



(c) 


(d) 


FIG. 2: The figure represents phase space trajectories, evolution of energy density and density parameter for 
interacting model II for the stable fixed point 3. The panel (a) corresponds to /3 = 0.7 and A = 1.3, for these 
values of the parameters we obtain = 0.76, Wtot = —0.57, = —0.75 and an accelerating attractor solution. 

The panel (b) corresponds to /3 = 0.3 and A = 1.3, correspondingly we get = 0.83, Wtot = —0.53, = —0.64 

and an accelerated attractor solution. In both the panels, black dots designate attractor stable point, and the 
stable point behaves as an attractive focus under the chosen parameters. The panels (c) and (d) are plotted the 
same values of the parameters as panel (a), and show scaling behaviours that gives an accelerated expansion. 


Above eigenvalues are negative for — oo < /3 < oo and — oo < A < oo (but A 0). 


The negativity of the eigenvalues exhibits the stability. This point is stable for all the values of /3 and 
A provided that A 7 ^ 0. We elaborate the autonomous system (IT71) numerically for the parameter choices 
(3 = 0.7, A = 1.3, and /3 = 0.3, A = 1.3, and obtain Wtot, W^ and as —0.57, —0.75, 0.76 and —0.53, 
—0.64, 0.83, respectively. The phase space trajectories of this stable point are shown in panels (a) and 
(b) of figure [U and the point behaves as an attractive focus. The lower panels of figure [5] exhibit the 
scaling behaviour that gives late time acceleration. The results of the stable point are abbreviated in 
table m 
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FIG. 3: The figure displays phase space trajectories, evolution of energy density and density parameter for 
interacting model III for the stable fixed point 3. The panel (a) corresponds to u = —0.3 and A = 1.57, for these 
values of the parameters we get = 0.65, Wtot = —0.46, = —0.71 and an accelerating attractor solution. 

The panel (b) corresponds to ct = 0.1 and A = 2.85, and correspondingly we obtain = 0.45, Wtot = 0.22, 
= 0.49, since total equation of state is positive therefore, attracting solution but not accelerating. In both 
the panels, black dots represent attractor stable point, and the stable point acts as an attractive focus under the 
chosen parameters. For panels (c) and (d), we use same values of the parameters as panel (a). Both panels are 
showing attractor behaviour that corresponds to scaling solutions. 


C. Interacting model III 

In this model we consider the coupling form as a linear combination of p'm and as 

A = a{pm + pV) 


( 20 ) 


For this interaction form we have 


H 

IP 

dX 

m 

dY 

dN 


3 

'1 + X^-Y^' 



2 

1 — a 



■iX + f-XY'‘ + 

1- 

'l + X^-Y"^' 

1 — a 




'1 + x^- 

1 — cr 


( 21 ) 

(22) 






























The total equation of state for this model can be obtained by using equations m and (EU as 


Wtot 


-1 + 


1 + X'^-Y 
1 — a 


21 




(23) 


This model has following stationary points: 


(1) X = 0, Y = 0 

In this case, the corresponding eigenvalues are 

Ml = - 2 i-i+a) < 0, for cr > 1 H 2 = 2 {-i%) < 0’ > 1/2 

This point is stable for cr > 1. 


( 2 ) X = VI - 2 cr, y = 0 

In this case, the corresponding eigenvalues are 

fj,i = 6 + < 0, for cr < 1/2 fj ,2 = 3 — y^ 3/2 — 3cr A < 0, for y^3/2 — 3cr A > 3 

This point is also stable for the above given conditions. 


(3) X = 


9-i-l+afY+rj 


~ 2\/3A 


2\/6A(3+(-l+cr)A2) ’ 

In this case, the corresponding eigenvalues are 


Ml = i (-3- 


12 


3r? 


24 r? 


3 ^ - 3(-l + cr)A^ + 

30 r) _I_ 18(7 7) 


V 2 ( 


9(-7+o-(2+13o-)) 

(- 1 +<t )2 


pi^ + 6aA2 + (-l + a)2A4-i7 


'(-l+<r)A2 (_l+cr)(3+(_l+cr)A2) (_l+cr)(3+(-l+<T)A2) 


,, _ 1 / o I 12 ‘T.l' 1 I I I ^“1 /9(-7+o-(2+13cr)) 

- 8 “ '1V--L + crjA + + y/Z - {=T+^ - 

24 ?7 30 77 I 18fT 77 

(-l+cr)A^ (-l+dr)(3+(-l+cr)A^) (-1+cr) (3+( - 1+cr) A^ ) 


^^ + 6aA2 + (-l + a)2A4 


where 

77 = y^(3 + (-1 + cr)A2)^ (9 - 6(1 + cr)A2 + (-1 + cr)2A'‘) 


(24) 


Above eigenvalues are negative for cr < 0.2 provided that cr V — 1 and A < 


6(l+cr) 

(- 1 +^) 2 - 


This point has both the parameters cr, A and shows the stability for the choices of the parameters 
cr < 0.2 (cr ^ — 1) and A < We numerically elaborate the autonomous system for the 

parameter choices cr = —0.3, A = 1.57, and a = 0.1, A = 2.85, and obtain Wtot, W^j), as —0.46, —0.71, 
0.65 and 0.22, 0.49, 0.45, respectively. We do not find accelerating solution corresponding to cr = 0.1 as 
it has positive equation of state. Also, we check numerically for all positive values of cr (0 < cr < 0 . 2 ) 
and do not find accelerating phase. The results of this stable point are concised in table HI The phase 
portraits of the point are displayed in panels (a) and (b) of figure [3l and the point acts as an attractive 
focus. The lower panels of figure [3] exhibit the scaling behaviour that provides late time acceleration. 
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TABLE I: We present the stable points 5, 3 and 3 for interacting models I, II and III respectively, these stable 
points have two parameters, namely, a, A and P, A and ct, A. We also introduce the expressions of 11,/,, Wtot and 
the condition of acceleration in terms of the parameters. Two numerical choices of the parameters are shown for 
each model. The symbols 5 and r] are used for shorting the expressions and given by equations m and iMll 
respectively. 


Model 

A 

Y 

Stable for 


Wtot 

IT,/, = 

Wtot 

Acceleration 

I 

^572 

y/3/2-3a 

a < -7/2, 

A > 

3 

OL 


for a < —1/2 

A{l-a) 

A(c-l) 

(1 —a) 

1 — OL 





a = -3.6 

A = 0.94 

0.73 

-0.78 

-1.06 

Yes 




a = —4.6 

A = 0.94 

0.60 

-0.82 

-1.36 

Yes 

II 

1 

^V[(1-/3^)(6 + 9/A2) 
-{1 + + S{1 + p 

-3/A= + 3P/X^)] 

for all values of 

P and A provided 
that A 7 ^ 0 

(3+a2)(1+/3)-5 

A^-3+.8(3+A^)-5 


A"+d(3 + A=) 
-5< 1 

2 V 6 A 

[(3 + A^) 

(l + /3)-5] 

2A2 

6 





P = 0.7 

A = 1.3 

0.76 

- 0.57 

-0.75 

Yes 




P = 0.3 

A = 1.3 

0.83 

- 0.53 

-0.64 

Yes 

III 

9-(1-<t) = A*+7, 

\J 9+6A2 — (1 —{T)2A4-j_7j 

cr < 0.2 7 ^ -1 

A< 

- 1 —<T 

9-(1-ct)^A*+7) 

1 


27-3(l-a)A= 

[4-(l-a)A2] 
—3 j/ < 0 

2\/^A(3-(l-a-)A2j 

2\f3\ 

< 

b 

1 

1 

18-6(1-<t)A2 

[(l-ajA^ 
(6-(l-a)A^) 
+r) - 9] 





cr = -0.3 

A = 1.57 

0.65 

-0.46 

-0.71 

Yes 




(7 = 0.1 

A = 2.85 

0.45 

0.22 

0.49 

No 


IV. CONCLUSIONS 

We studied interaction of quintessence with dark matter in spatially flat Universe. In the absence 
of fundamental theory of specific interaction in the dark sector, the choice of interaction strength in 
the conservation of energy equations was phenomenological and heuristic. In this paper we considered 
three phenomenological interacting quintessence cosmological models as ap'm, Pp'rp, and a{p'm + P 4 ,)- 
Our primary object was to inspect whether there exist late time accelerated scaling attractor having 
^de/^dm of the order one. We studied dynamical behaviour and phase space analysis of the models 
under consideration. We focussed on the stable points which could give rise to scaling attractors. In 
all the models we obtained fundamental cosmological observables like Wtot, W^ corresponding to 
two numerical choices of the parameters. For the interacting model I we found that the fixed point 5 
is stable for a < —7/2 and A > y^2/3. The phase space trajectories, evolution of energy density and 
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density parameter for different numerical choices of the parameters are shown in figure [TJ The fixed 
point 5, in this case, corresponds to accelerated scaling attractor with ^Ide/^dm = 0(1)- The point 3 
of interacting model II shows stability for all values of (3 and A provided that A 0. The phase portraits, 
evolution of energy density and density parameter for different numerical values of the parameters are 
displayed in figure[5] Clearly, this is scaling solution with the required property. For interacting model III 
we noticed that the fixed point 3 is stable for tr < 0.2 —1 and A < y^6(1 + cr)/(I — cr). Figure 12] shows 

the phase space trajectories, evolution of energy density and density parameter for different values of the 
parameters. Our analysis shows that accelerating attractor, in this case, exists only for negative values 
of (j. The lower panels of figure [2] shows accelerated attractor solution for cr = —0.3 and A = 1.57. For 
all the models, we obtained late time accelerated scaling attractor having Ude/^dm = 0(1). Therefore 
all the models considered in this paper are viable to solve the coincidence problem. 
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